NOTE ON RADIAL DUNKL PROCESSES 



NIZAR DEMNlQ 



Abstract. This note encloses relatively short proofs of the following known 
results: the radial Dunkl process associated with a reduced system and a strictly 
positive multiplicity function is the unique strong solution for all time t of a 
stochastic differential equation of a singular drift (see [T2] for the original proof 
and [1] for a proof under additional restrictions), the first hitting time of the 
Weyl chamber by a radial Dunkl process is finite almost surely for small values of 
the multiplicity function. Our proof of the second mentioned result gives more 
information than the original one. 

1. Preliminaries 

To be self-contained, we begin by pointing out some facts on root systems and 
radial Dunkl processes. We refer to [TT] for the Dunkl theory, to [8] for a back- 
ground on root systems and pE] and references therein for facts on radial Dunkl 
processes. Let (V, (, )) be a finite real Euclidean space of dimension m. A reduced 
root system R is a finite set of non zero vectors in V such that : 

1 R n Ra = {a, -a} for all a G R, 

2 a a (R) = R, 

where a a is the reflection with respect to the hyperplane H a orthogonal to a: 



a a (x) = x — 2 x ' ' a, \a\ 2 := (a, a) x £ V. 



a 



A simple system S is a basis of span(i?) which induces a total ordering in R. A 
root a is positive if it is a positive linear combination of elements of S. The set 
of positive roots is called a positive system and is denoted by R + . The (finite) 
reflection group W is the group generated by all the reflections a a for a G R. 
Given a root system R with positive and simple systems R + , S, define the positive 
Weyl chamber C by: 

C := {x G V, (a, x) > OVa G R + } = {x G V, (a, x) > OVa E S 1 } 

and C, dC its closure and boundary respectively. The radial Dunkl process X is 
defined as the C-valued continuous paths Markov process whose generator is given 
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by: 



^ k u(x) = -Am (x) + V k(a)^-. Y 1 



where u G C^(C) satisfying the boundary conditions (V«(i),«) = for all x G 
H a , a G R + , and k(a) > is a multiplicity function (a H^-invariant function). 

2. Motivation 

In order to motivate the reader and prepare for the first result, we exhibit some 
known examples. The first one is given by the rank-one case. It corresponds to 
V = R, R = B 1 = {±1}. Hence k(a) := k > and X is a Bessel process ([10]) of 
index v — k — 1/2. When ki > and X = x > 0, it is the unique strong solution 
of the following stochastic differential equation of singular drift: 

k 

dX t = dB t H dt, t > 0, 

Xt 

where B is a standard Brwonian motion. A multivariate well known example is 
given by the A m _i-type root system defined by: 

A m _i = {±{ei - ej), 1 < i < j < m}, 

with positive and simple systems given by : 

R + = {ti — ej, 1 < i < j < m}, S = {e« — 1 < i < m}, 

where (ei)i<j< m is the standard basis of M m . In this case, V = M m , the span of 
R is the hyperplane of lR m consisting of vectors whose coordinates sum to zero 
and C = {x G lR m ,a;i > ••• > x m }. Besides, there is only one orbit so that 
k(a) := k > and X = (X l )i<i< m satisfies : 

(1) dXl = dp] + k^ — — r Ki<m } t<r 

j^i t t 

with Xq > ■ • • > X™, where [y % )i are independent BMs and r, the first collision 
time, is defined 

t := inf{t, X\ = X{, for some (i, j)}. 

This process was deeply studied in [2] and it was shown that (TjQ) has a unique 
strong solution for all t > and Xq > • • - > X™ provided that k > 0. Then, it 
was proven in [12] that the radial Dunkl process is the unique strong solution of 



dX t = dB t - V$(X t )dt, X eC 

where $(x) = — J2 a eR+ k( a ) ln((a, x)) subject to k(a) > for all a G R. Schapira's 
proof relies on a similar existence and uniqueness result for radial Heckman-Opdam 
processes and an auxiliary result showing that radial Dunkl processes are limiting 
(in law) processes of rescaled radial Heckman-Opdam processes. Meanwhile and 



independently, Chybiryakov provides another more restrictive proof based on mar- 
tingale problems however assuming that the simple system is a basis of V which 
is not valid for R = A m _x, and that the starting point X G C. Both proofs are 
long and in |3], the result was stated in Remark 15 p. 13. Here we give a rela- 
tively short proof based on a result by Cepa and Lepingle on stochastic differential 
equations of singular drifts ([2]) and an additional lemma which roughly says that 
if the process reaches the boundary of the Weyl chamber, then the inward normal 
vector at that time can not be zero (it will instantaneously reflect). 

Next, we interest in the first hitting time T of dC and prove, via stochastic 
calculus, the fact that T < oo a.s. when < k(a) < 1/2 for at least one a G R + . 
This result was proved in [I] using local martingales and the proof is technical and 
long. Our proof provides more information and is a generalization of the one given 
in [2]. Indeed, for such an a, we prove that the process (a,X) hits zero a.s. by 
comparing it to a Bessel process of dimension in the interval [1, 2[. 



3. Radial Dunkl Process : Existence and Uniqueness of a strong 

solution 

Theorem 1. Let R be a reduced root system and recall that: 

— — k(a)\n((a,x}) : = V} k(a)0({a,x)), x E C, 

where k(a) > for all a G R + . Then the SDE 

(2) dX t = dB t - V<S>{X t )dt, Xq G C 

where X is an adapted continuous process valued in C and B is a Brownian motion 
in W 71 , has a unique strong solution. 

Proof. From Theorem 2. 2 in [2], we deduce that: 

(3) dX t = dB t - V$(X t )dt + n(X t )dL t , X eC 

where n(x) is a (unitary) inward normal vector to C at x , L is the boundary 
process satisfying: 

dL t = l{x t £dC}dL t , 
has a unique strong solution for all t > 0. Moreover: 



(4) 
(5) 



E 



E 



-{Xt&dC} 



dt 



\V$(X t )\dt 



0. 



< oo 



for all T > 0. Thus, it remains to prove that the boundary process vanishes. To 
proceed, we need two Lemmas. 
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Lemma 1. Set dG t '■— n(X t )dL t . Then, for all a G R+, 

l{(x t ,a)=o}(dG t7 a) = 0. 

Proof. The proof is roughly an extension to arbitrary root systems of the one 
given in [2J for R = An-i . In order to convince the reader, we provide an outline. 
On the one hand, the occupation density formula yields: 

rt 

m(a,X))9 (a)da= (a, a) / d'((a,X s ))ds 

Jo 

where L^((a,X)) is the local time up to time t at the level a > of the real 
continuous semimartingale (a, X) > ([1Q]). On the other hand, the following 
inequaliy holds (instead of (2.5) in [3]) for all a G C: 



(V$(x), x — a) = k{a)9 ((a, x))(a, x — a) 



(i) 



> k(a)[b a 9 ((a, x)) — c a (a, x — a) — d a ] 



> min (b a k(a)) 9 ((a, x)) — \x — a\ k(a)c a \a\ — k(a)da 

:= A 9'((a,x)) - B\x - a\ - C 

by Cauchy-Schwarz inequality, where in (1), we used equation (2.1) in [3j: let g 
be a convex C 1 -function on an open convex set D C M m , then Va G D, there exist 
b, c, d > such that for all x G D: 

(Vg(x),x — a) > b\Vg(x)\ — c\x — a\ — d. 

Note also that A > since b a k(a) > for all a G R+. Then, the continuity of X, 
([5]) and the above inequality yield: 



/ 9' ((a,X s ))ds < oo 
Jo 



which implies that: 

POO 

m(a,X))9'(a)da < oo 



Thus, Lt((a, X)) = since the function a \—> 9' (a) is not integrable at 0. The next 
step consists in using Tanaka formula to compute 

dZ t :=d[(a,X t )-((a,X t )) + ] 

= l{{ a ,xt)=o}(oc, dB t ) - l{(a,x t )=o}(a, V$(X t ))dt + l{^ Xt )=o}(a,dG t ) 
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for a G A. It is obvious that the second term vanishes. The first vanishes too since 
it is a continuous local martingale with null bracket (occupation density formula). 
As X t G C, then dZ t = a.s. which gives the result. ■ 

Lemma 2. Let x G dC. Then (n(x), a) 7^ for some a 6 5 such that (x, a) = 0. 

Proof: Let us assume that (n(x), a) = for all a G S such that (x, a) = 0. The 
idea is to find a strictly positive constant e such that x — en(x) G C then use the 
definition of inward normal vectors (see below) to conclude that n(x) = 0. Our 
assumption implies that (x, a) > for all a G S such that (n(x),a) 7^ 0. If such 
simple roots do not exist, that is, (n(x), a) = for all a G S, then x — en(x) G C 
for all e > 0. Otherwise, if (n(x), a) < for these simple roots, then x — en(x) G C 
for all e > 0. Finally, if none of these conditions is satisfied, choose 

(x, a) 

(J < e < mm — -, 

<x,a)>0,{n.(a:),o)>0 (n(x),a) 

to see that x — en(x) G C Substituting a = x — en(x) in the definition of the 
inward normal n(x) to C at x 

(6) (x - a,n{x)) < 0, a G C, 

it follows that n(x) is the null vector which is not possible. ■ 
Now we proceed to end the proof of Theorem [H If X t G H a for one and only one 
a G A. Then, n(X t ) = a/\a\ and Lemma [1] gives 

l{(Xi,a)=0}(^G' t , a) = l{{x t ,a)=0}\oc\dL t = 0. 

Hence, (L t ) t >o vanishes a.s.. More generally, Lemma [2] asserts that 

{X t G dC} C U QeA {(X t , a) = 0, (n(X t ), a) ^ 0} 
for all t. It follows that 



< L t < [ 



l{{X t ,Q>=0,(n(X t ),a)^0}0?-^s 
1 



l{{X t ,a)=0,{n(X t ),a)^0}(^(^s), Ol)dL s 



(n(X s ),a) 



by Lemma [U 



Remarks. 1/When m = 1, (JT t ) t > is a Bessel process of dimension 5 = 2k + 1 
and ko > 5 > 1. It is well known that the local time vanishes (see Ch. XI in 

HD]). 

2/Since the boundary dC of C is a cone, then for x G dC, one has ex G dC for all 
c > 0. Choosing ai = Cix for < Ci < 1 and 02 = C2X for C2 > 1, one gets from 
the definition of normal vectors that (n(x), x) = so that n(x) G x -1 . 
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3.1. Another example and an open question: the £> m -type root system 
and /3-Laguerre processes. The _B m -type root system is defined by 

R = {±6j, 1 < i < m, ±6j ± ej, 1 < i < j < m}. 

Its positive and simple systems are given by 

R + = {ei, 1 < i < m, ei±ej, 1 < i < j < m}, S = {e^ — e^+i, 1 < i < m — 1, e m }. 

It is known that the radial Dunkl process of type A fits the eigenvalues process 
of the symmetric and the Hermitian Brownian motion (|7J). The B-type root 
system turns out to be related to the eigenvalues process of Wishart and Laguerre 
processes (see [I], [5]) which is the unique strong solution of: 

\(t) + X k (t) 



(7) d\i(t) = 2^X~^)dv i (t)+(3 



k^i 



' \(t) - X k (t) 



dt, 1 < i < m 



for (5 — 1, 2 and S > m + 1, m respectively, where {vj)i are independent BMs and 
Ai(0) > • • • > A m (0) > 0. The range of time is given by t < r A Ro where r is as 
before and 

R := inf{t, X m (t) = 0} 

is the first hitting time of zero. It is known that r = oo almost surely and that 
Ro = oo almost surely if 8 > m + 1, m respectively (p],[5]). 

Define the /?-Laguerre process as a solution, whenever it exists, of 
arbitrary j3, 5 > up to R A r. Set := \/A7, then, for t < t A Rq: 



with 



driit) = dviit) + 



dvi(t) + 



2ri(t) 
(3(5 - 



rf + r| 



dt 



m + 1) 



2n(t) 



1 , ^ 



dvAt) H j-rdt + ki 

nit) z — ' 



1 



n(t) -r.it) 

1 



+ 



+ 



Vn{t)-r 3 {t) n(t)+rj(t) 



ri(t)+rj{t) 



dt 



dt 



with 2k = (3(5 — m + 1) — 1, 2ki = (3. Consequently, the process r = (r 1? . . . , r m ) 
is a S m -radial Dunkl process and using Theorem [TJ one claim that (r(t)) t > the 
unique strong solution for all t > of ([7]) provided that k , k\ > 0, that is, 
(3 > 0, 5 > m — 1 + (1/(3). This improves results from matrix theory: for Wishart 
processes, take (3 — 1 so that the strong uniqueness for the eigenvalues process 
holds for all 5 > m. For Laguerre processes, take (3 = 2 so that strong uniqueness 
holds for 5 > m — 1/2. Actually, the range of dimensions 5 for which the strong 
uniqueness holds is extended for the eigenvalues process and one can even start at 
t = from A(0) e C (see Theorem 2 p. 561 in [5]). However, it is known that the 
eigenvalues process of a Laguerre process is the unique solution in law of ([7]) when 
5 > m — 1 provided that Ai(0) > • ■ ■ > A m (0) (Theorem 4 p. 562) and a similar 
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result holds in the real case (PQ). Besides m — 1 < 5 < m, (3 = 2 <=>■, —1/2 < k < 
0,ki = 1 and it is already known that Bessel processes (rank one case) are well 
defined for < 5 < 1 ^ — 1/2 < k < via excursion theory (|13j). Thus, one 
wonders if the radial Dunkl process can be defined as the unique solution in law of 
(j2J) for non positive k (note that k = —1/4,-1/2 are singular values for the Dunkl 
intertwining operator, [6]). 



4. FlNITENESS OF THE FIRST HITTING TIME OF THE WEYL CHAMBER 

Let T := inf {t > 0, X t G dC} be the first hitting time of the Weyl chamber. In 
[2], where R = An-i and T = r = inf{t > 0, X\ = X{ for some authors 
showed that To < oo a.s. if < k± < 1/2. More generally, it was shown that if 
< k(a) < 1/2 for some a G S, then To < oo a.s. Here we prove a more precise 
statement: 

Proposition 1. Let oq G S and T ao := inf{t > 0, (ao,X t ) = 0} such that To = 
m fa eAT ao . If < k(a ) < 1/2, then ((a ,X t )) t >Q hits a.s. 0. Therefore T < 
T ao < oo a.s. 

Proof. Assume k(a) > for all a G R and let a G S. Our scheme is a 
generalization of the one used in [2], thus we shall show that the process < a , X > 
is a.s. less than or equal to a Bessel process of dimension 2k(a>o) + 1 := 2&o + 1. 
The result follows from the fact that 2k(ao) + 1 < 2 when k(a) < 1/2 so that the 
Bessel process hits zero a.s.. Using (|2]), one has for all t > 

d(a ,X t ) = 1^01^74+ y~] k(a) | Q ' dt 

= \a \d^ t + h dt + S~] k(a) , a ' "°. dt. 

{at ,X t ) (a, At) 

aeR+\a 



Set 



R = U p j=1 R j , 



where Ri, 1 < j <p denote the conjugacy classes of R under the H^-action, then 

r + = uUK 

so that: 



I 12 P I \ 

, «o (a.ao) 
d<a ,X t >=\a \dj t + k - —rdt+ykj } j -. — z^rdt. 



For a conjugacy class Ri and a G R? , if (a, a ) := a(a) > then, it is easy to check 
that (o"o(a), cto) = —a(a) where cr is the reflection with respect to the orthogonal 
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hyperplane H ao . Note that cr Q (a) belongs to the same conjugacy class of a and 
that a (a) G R+ \ a o f° r a G R+ \ a o ( see Proposition 1. 4 in [S]). Hence, 

a(a , -Xt) = «o d 7t + «(n ttt^ - > % > - — dt. 

(a ,JQ) fr? V (a,^t) (o"o(a),Xt) 



a(a)>0 



Furthermore, 



(a,a ) . , . . , , (a ,X t ) 

a — cr (aj = 2- r-a =^ (a — o- {a), X t ) = 2a(a) 



(a ,a ) ' 1 ' |a | 



Consequently: 



12 



ao 

d(ao) -X't) — |ao|^7t + k - Trvdt + F t dt 



(oto,X t ) 

where F t < on {T ao = oo}. Using the comparison Theorem in [9] (Proposition 
2. 18. p. 293 and Exercice 2. 19. p. 294), one claims that (a ,X t ) < Y^ H 
for all t > on {T ao = oo}, where Y x is a Bessel process defined on the same 
probability space with respect to the same Brownian motion, of dimension 2k + 1 
and starting at Y = x > (ao,X ) > 0. This is not possible since a Bessel process 
of dimension < 2 hits a.s. ([10]. Chap. XI). ■ 
Open question: Given two simple roots a,\,a.2 G S such that 

< fc(ai) ^ fe(aa) < 1/2, 

that is belonging to different orbits, we already know that T a . < oo,i — 1, 2 a.s.. 
Can we compare T ai to T Q2 ? This may be possible if one finds two processes Ri, R 2 
such that for instance 

F(( ai ,X t ) < R^t) < (a 2 ,X t ) < R 2 (t), for all t > 0) = 1 

and i?2 hits zero a.s. 
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